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ABSTRACT 

We present a simple physical mechanism that can account for the observed stellar 
mass spectrum for masses M* > 0.5Mq. The model depends solely on the competitive 
accretion that occurs in stellar clusters where each star's accretion rate depends on 
the local gas density and the square of the accretion radius. In a stellar cluster, there 
are two different regimes depending on whether the gas or the stars dominate the 
gravitational potential. When the cluster is dominated by cold gas, the accretion 
radius is given by a tidal-lobe radius. This occurs as the cluster collapses towards a 
p cx R~ 2 distribution. Accretion in this regime results in a mass spectrum with an 
asymptotic limit of 7 = —3/2 (where Salpeter is 7 = —2.35). Once the stars dominate 
the potential and are virialised, which occurs first in the cluster core, the accretion 
radius is the Bondi-Hoyle radius. The resultant mass spectrum has an asymptotic 
limit of 7 = — 2 with slightly steeper slopes (7 w —2.5) if the stars are already mass- 
segregated. Simulations of accretion onto clusters containing 1000 stars show that 
as expected, the low-mass stars accumulate the majority of their masses during the 
gas dominated phase whereas the high-mass stars accumulate the majority of their 
massed during the stellar dominated phase. This results in a mass spectrum with a 
relatively shallow 7 « 3/2 power-law for low-mass stars and a steeper, power- law for 
high- mass stars — 2.5 < 7 < —2. This competitive accretion model also results in a 
mass segregated cluster. 



Key words: stars: formation - stars: luminosity function, mass function 
clusters and associations: general. 
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1 INTRODUCTION 

One of the most fundamental unsolved problems in astron- 
omy is the origin of the distribution of stellar masses, the so- 
called Initial Mass Function (IMF) . Significant observational 
work has gone into establishing the exact form of the mass 
spectrum (c.f. Scalo 1986; Kroupa, Tout & Gilmore 1990; 
Kennicut 1998) and investigating it's universality (Kenni- 
cut 1998; Leitherer 1998; Scalo 1998). Unfortunately this 
work has not led to an increase in our understanding of 
how the mass spectrum originates. This shortcoming is in- 
trinsically linked to our lack of understanding of the star 
formation process. Thus, theories of the IMF have generally 
involved sufficient free-parameters to ensure that they can 
reproduce the general form of the mass spectrum but also, 
unfortunately, ensuring that they do not help us discern the 
physics involved. 

Observations of star forming regions have clearly shown 
that most stars are formed in clusters (e.g. Lada et al. 1991; 
Lada, Strom and Myers 1993; Meyer & Lada 1999; Clarke, 
Bonnell & Hillenbrand 2000). Furthermore, the individual 
stellar masses found in these clusters span the observed 



range of stellar masses and are consistent with being drawn 
from the field star IMF (Hillenbrand 1997; Hillenbrand & 
Carpenter 2000). Thus, the basic physics which results in 
the observed mass spectrum is present in this fundamen- 
tal building block of the star formation process: the young 
stellar cluster. Competitive accretion in gas-rich stellar clus- 
ters results in a spectrum of stellar masses (Bonnell et 
al. 1997, 2000) and in this paper we investigate whether 
this process can lead to the observed IMF. 

There have been several theories advanced to explain 
the observed distribution of stellar masses (cf. Clarke 1998; 
Meyer et al. 2000). These theories generally involve either 
fragmentation or accretion plus feedback as the basic physics 
in setting the stellar masses. Fragmentation of a giant molec- 
ular cloud can set the stellar mass distribution if the CMC 
clump-mass distribution is, for some reason, of the same 
form as the IMF (Elmegreen 1997), or if feedback changes 
the basic fragment (Jeans) mass (Silk 1977). Alternatively, 
accretion onto a small isolated protostellar core can set the 
final masses but this requires a halting mechanism (eg feed- 
back) that really determines the mass distribution (Adams 
& Fattuzzo 1996). The main criticisms of these theories are 
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that a) they rely on additional processes (e.g. feedback or an 
initial mass distribution) and b) that they neglect the fact 
that the stars are not formed in isolation but in clusters and 
therefore can interact during their formation. Furthermore, 
they are generally unable to explain why clusters are ini- 
tially mass-segregated (eg, Hillenbrand & Hartmann 1998; 
Bonnell & Davies 1998) which should be a requirement of 
any theory that attempts to explain the IMF. In this paper, 
we discuss how competitive accretion in clusters can set the 
distribution of stellar masses and result in a mass-segregated 
cluster. 



2 ACCRETION IN CLUSTERS 

Young stellar clusters are generally found to contain signifi- 
cant amounts of gas. In fact, the total gas mass is typically 
many times the mass in stars (Lada 1991). This gas is a rem- 
nant of the star formation process and can be understood 
in terms of the inefficiency of fragmentation where the ma- 
jority of a system's mass is not immediately included in the 
fragments but is accreted by them over longer time periods 
(Larson 1978; Boss 1986; Bonnell & Bastien 1992; Burkert 
& Bodenheimer 1993; Klessen, Burkert & Bate 1998). This 
occurs due to the non- homologous nature of gravitational 
collapse (Larson 1969) so that the amount of mass that has 
reached the region where and when fragmentation occurs, is 
small. Thus, the initial fragment mass should generally be 
small, M « Mq. Furthermore, in a cluster environment, 
the fragment mass, if it is proportional to the Jeans mass, 
should be smallest in the centre of the cluster where the den- 
sity is highest, in direct contradiction to the observed mass 
segregation (Bonnell & Davies 1998; Bonnell et al. 1998). It 
is thus unlikely that the mass distribution is solely deter- 
mined by the fragmentation process. 

Accretion of the residual gas in clusters can play a cru- 
cial role in setting the final stellar masses (Zinnecker 1982; 
Larson 1992). As all the stars are embedded in the same en- 
velope of gas, they compete for the gas as they move through 
the cluster. This competition results in unequal mass accre- 
tion rates and hence a mass spectrum 

Simulations of accretion in small stellar clusters (Bon- 
nell et al. 1997) have shown that even if the stars have equal 
initial masses, they accrete unevenly with those stars near 
the centre accreting more than those near the outer parts of 
the cluster. This dependency of the accretion rates on posi- 
tion in the cluster is even clearer in larger clusters (N = 100; 
Bonnell et al. 2000) where accretion rates can vary by over 
an order of magnitude from the outside to the cluster centre 
when the cluster is centrally-condensed. The physics behind 
this differential accretion is the effect of the cluster potential 
in funneling material down to the centre plus the individ- 
ual star's ability to capture this material. Initially uniform 
clusters also accrete differentially as they evolve towards a 
centrally-condensed configuration, due to the local compe- 
tition between stars. The final cluster displays a significant 
degree of mass segregation due to the differential accretion 
(Bonnell et al. 2000). 

A star's accretion rate depends on it's cross section for 
accreting material, nRl cc , the gas density through which it 
moves, p and the velocity with which it moves relative to 
the gas, Vr e v. 



M* = TTptW&e. (1) 

The question is what is the accretion radius, -Race- An iso- 
lated star of mass M*, moving in a uniform medium with 
velocity V Te \ relative to the gas with sound speed c s , has an 
accretion radius given by Bondi-Hoyle accretion where 

Rbh = 2GM«/(V^ei + c?s) (2) 

(e.g. Bondi & Hoyle 1944; Ruffert 1996). On the other hand, 
in a cluster environment where cold gas dominates the po- 
tential, the accretion radius can be described by the star's 
tidal-lobe radius. 

where M* is the star's mass, M enc is the cluster mass interior 
to the star's position in the cluster R (e.g. Paczynski 1971; 
Frank, King & Raine 1985). Material inside this tidal- lobe 
is more attracted to the star whereas outside this lobe the 
acceleration is dominated by the cluster potential. 

Simulations of accretion in gas-dominated stellar clus- 
ters have shown that the accretion is well modelled by a 
tidal-lobe radius as long as the gas dominates the gravita- 
tional potential (Bonnell et al. 2000). When this is the case 
and the gas is unsupported, as expected in order for the clus- 
ter to have fragmented, both the gas and the stars collapse 
towards the cluster centre. The stars are impeded from viri- 
alising due to the changing potential, and the gas accretion. 
Therefore the gas and stellar velocities are locally correlated 
such that the relative gas velocity is small. This results in a 
large -Rbh such that Rbh > Rtidai and therefore the smaller 
Rtidai is the appropriate accretion radius. 

In contrast, once the stars dominate the gravitational 
potential, which occurs first in the core of the cluster due 
to the higher accretion rates there and the sinking of the 
more massive stars, the stars virialise and their velocities 
relative to the gas become large. When this occurs the Bondi 
Hoyle radius becomes much smaller due to the large gas 
velocities such that Rbh < -Rtidai- At this point, the Bondi- 
Hoyle radius is a better description of the accretion radius 
(Bonnell et al. 2000). 



3 THE INITIAL MASS FUNCTION 

In order to use the above formulations of the accretion rate 
to determine the resultant mass spectrum, we need a basic 
model of a stellar cluster. The cluster is initially gas dom- 
inated in order to provide the mass reservoir for accretion. 
Observations of the youngest clusters estimate that up to 90 
per cent of the mass is in gas (Lada 1991). Such systems are 
unlikely to be in equilibrium, and regardless of their initial 
configurations, will collapse towards a 

p oc RT 2 (4) 

density distribution (Larson 1969; Hunter 1977). The stars 
will also collapse, even if initially virialised, due to the 
change in the potential. Thus we assume that the stars fol- 
low the same distribution: 

n oc R~ 2 (5) 

This profile is a good approximation to the stellar distribu- 
tion of most clusters (Binney & Tremaine 1987) including 
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the youngest clusters that have been studied such as the 
Orion Nebula cluster (Hillenbrand & Hartmann 1998). Fur- 
thermore, it is the expected distribution produced by a vi- 
olent relaxation of a non-equilibrium purely stellar cluster 
(Lynden-Bell 1967). 

The gas distribution will be modified by the accre- 
tion and hence removal of material. The accretion rates are 
largest in the centre and hence it is there that they will first 
affect the density profile. This gas removal destabilises the 
equilibrium of the isothermal sphere allowing gas to infall 
from further out in the cluster. In the case of infall onto 
a point-mass, the density profile then takes on the form 
(Hunter 1977; Shu 1977; Foster & Chevalier 1993) 



p oc R 2 . 



(6) 



We thus assume an initial density profile of an isothermal 
sphere (equation ^) which evolves towards a profile com- 
posed of the shallower, equation^ profile in the central 
regions while the outside maintains the isothermal sphere 
profile. The transition between these two regimes occurs 
when gas accretion has increased the stellar masses suffi- 
ciently that the stars dominate the gravitational potentialin 
the cluster centre. We discuss the asymptotic limits to the 
accretion and mass spectrum in these two cases below. 



3.1 Gas dominated potentials 

In addition to the gas density we require estimates of the 
relative gas velocity and of the accretion radius. Where gas 
dominates the potential (p oc R~ 2 ), the accretion radius is 
given by the tidal-lobe radius (Bonnell et al. 2000) 



Rn 



Rtid 



(7) 



In order to evaluate the tidal-lobe radius we need the mass 
enclosed at a particular radius in the cluster, R, 

R 

(8) 



M cnc = / 4irpr 2 dr, 
Jo 

which in a p cx R~ 2 distribution is 

M enc tx R. 



(9) 



Simulations of accretion in stellar clusters have shown 
that in gas-dominated clusters, the relative gas velocity is 
generally subsonic as both the stellar and gas velocities are 
dominated by the acceleration of the cluster potential (Bon- 
nell et al. 2000). Alternatively, if the stars are virialised then 
the relative gas velocity is the velocity dispersion, 



R 



(10) 



Thus, if the gas distribution follows equation(^), then we 
have 



V„ 



constant. 



(11) 



In either case the relative gas velocity is independent of po- 
sition in the cluster. 

We can now use the above to consider how accretion in 
a stellar cluster results in a mass spectrum of the form 



dN oc M?dM* 



(12) 



where 7 = —2.35 is the Salpeter IMF. First of all, we note 
that from equation(^) the number of stars at a given radius 
is constant: 

dN = n{R)A-KR 2 dR cx const x dR. (13) 
Starting from the general equation for the accretion rates, 



M, = ■npV tA Ri 



(14) 



and adapting it for the case where gas dominates the poten- 
tial such that p oc R~ 2 and i? aC c ~ -Rtidai, then we have 

2 



oc g{t)R~ 2 x 



• M, 1 3 
.Af- 



(15) 



The term g(t) includes all time-dependence of the homolo- 
gously evolving cluster and R is the initial stellar position 
in the cluster. 

Integrating equation(^B|) with respect to time, t, and in 
the limit where accretion dominates the final mass, we get 

M» oc R- 2 h(t) 3 , (16) 

where h(t) — J g(t)dt. Thus, at a time ianal, when the ac- 
cretion is halted, we have 

RotM'^ (17) 
and thus 

dRcx M'^dM,. (18) 
This results in a mass spectrum (equation |l3| ) 

dN oc M~^dM t . (19) 

In general, the p oc R~ 2 profile will not extend right to 
the centre of the cluster. Instead, the central regions must 
have a flatter profile to stop the density from becoming infi- 
nite. A collapsing isothermal gas cloud has a p oc R~ 2 den- 
sity profile in the outer regions and an ever decreasing cen- 
tral region with near uniform density (Larson 1969). Stars 
in this region will see a shallower density profile which will 
affect their accretion rates (equation hj). Adapting equa- 
tion([l5]) for a density profile p oc R~ a while maintaining a 
uniform V te \ gives 



M» oc R~ a 

and analogous to equation (H), 



dR oc M» 



dM,. 



(20) 



(21) 



Taking the stellar density to be proportional to the gas den- 
sity, n(R) oc p(R), equation (p^) becomes 



dN 



oc M, 
oc M, 



n{R)AirR dR oc R 2 ~ a dR, 

-2. -±_1 



dM* 



(22) 



"dM,, 



and 7 is now 7 = — 3/a. From equation(K3) we see that 
when the density profile is shallower than R~ , the resulting 
mass spectrum is steeper than M^ 3 ^ 2 . Figure |l| plots the 
resulting IMF using equation (^) and a density profile that 
combines an outer region with a p oc R~ 2 profile with an 
inner Gaussian density profile (fitted onto the outer region 
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Mass 

Figure 1. The resultant IMF as derived from equation(p^), 
is plotted in log mass bins for a gas dominated cluster with 
p oc R~ 2 on the outside and a Gaussian central density profile 
(lower panel). The masses are scaled such that the average mass 
is OAMq. The upper panel shows the corresponding value of 7 
where dN = M^dMt. The lower-mass stars have a 7 = —3/2 
(r = —1/2) mass spectrum (where Salpeter is 7 = —2.35) due to 
the p oc R~ 2 whereas the more massive stars have an ever steep- 
ening mass spectrum due to the central Gaussian density profile 
(upper panel). 



where it has the same slope). The mass spectrum is plotted 
as a function of the logarithm of the mass such that the slope 
is r = 7+I. The low-mass stars that form in the outer region 
have the 7 = —3/2 mass spectrum whereas the higher-mass 
stars formed in the inner regions have an ever steeper mass 
spectrum due to the central Gaussian density profile that 
converges to uniform density. It is also worth noting that 
if the collapse gives a density profile that is steeper than 
R~ 2 as was found in the outer regions of the collapse of 
a Bonnor-Ebert spheroid (Foster & Chevalier 1993), then 
this translates (equationf22]) into a shallower IMF for the 
low-mass stars that are in this region. 



3.2 Stellar dominated potentials 

In contrast to a gas dominated system, when the cluster 
potential is dominated by the stars and they are virialised, 
the appropriate accretion radius is the Bondi-Hoyle radius 
(Bonnell et al. 2000). The Bondi-Hoyle radius models better 
the accretion radius as the gas and stellar velocities are no 
longer correlated. The gas is primarily infalling from further 
out in the cluster onto the stellar dominated core. When 
infall is onto a point mass, the resulting velocity profile is 



Kcl oc R 



-1/2 



(23) 



which assuming a constant mass infall rate with radius im- 
plies a density profile of the form of equation (^), 



p oc R 



-3/2 



(24) 



In the case of accretion onto a stellar cluster, the distributed 
mass and the presence of multiple accretion sources will 
modify these profiles. We can generalise them as 



Kci oc R~ 
and 

p oc -R -< °, 



(25) 



(26) 

under the constraint that the mass infall, AnR 2 pV rc \ de- 
creases with decreasing radius, 77 + £ < 2. we have 



oc V le \pR% n 



oc p- 



(27) 



oc g(t)R- ( M 2 R 3v 
oc g{t)M*R 3n -S. 



As above, g(t) includes the time-dependence of the homolo- 
gously evolving cluster and R is the initial stellar position in 
the cluster. In the case where 3r/ — £ = as occurs for accre- 
tion onto a point mass, the accretion rates are independent 
of position in the cluster. 

Assuming that the cluster maintains the same form 
(377 — C ~ constant), solutions to equation (p^) are 

M»(0) 



M, 



-/3M«(0)h(t) 



(28) 



where M*(0) is the initial stellar mass, h(t) — J g(t)dt, and 
j3 is basically the accretion clock which depends on the gas 
content and on R 3rl ~ ( but not on time. The primary differ- 
ence between this case and the preceding tidal-lobe accretion 
is the higher dependency on the stellar mass such that the 
initial mass is always important. This is the case for any 
accretion rate which depends on the stellar mass to a power 
greater than one. 

The resulting mass spectrum from such an accretion 
rate was studied by Zinnecker (1982) and we follow his 
derivation here. As there is a direct correlation between ini- 
tial and final mass, then there is a direct mapping from the 
initial mass spectrum to the final one, 



F(M*)dM, = F(M»(0))dM»(0). 



(29) 



Inverting equation (]2S|) and inserting it into equation (|29j), as 
long as f3 is independent of the initial masses, we have 



F(M») 



nM.(o))(^)" 

( M *(0)) (^Toy- 



M. (0) 



PM,(0)h(t) 
2 



(30) 



F(M.(0))(l-/3M.(0)ft(t))', 



= F 



1+0M, 



h(i)) ( 



M, (0) 1 
M. ) 



In the limit (3M t h(t) >> 1 (M» >> M»(0)) this gives 
F(M t ) oc M~ 2 , 



(31) 



and thus as long as there is an initial (small) range of stellar 
masses, the asymptotic limit of the mass spectrum is 



dN oc M7 2 dM t 



(32) 



This is somewhat complicated if there is a correlation 
between the initial masses and radii in the cluster. In this 
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Figure 2. Simulated IMFs for a combined 28 clusters contain- 
ing 1000 stars each are plotted in log mass bins. The clusters 
are assumed to have stellar dominated potentials and accrete ac- 
cording to either equation(^vj) with no spatial dependence (up- 
per panel) or equation(|33[) for an initially mass-segregated cluster 
(lower panel). The stars are in a n oc R~ 2 configuration and are 
initially a delta function in mass (upper panel) or the expected 
7 = —3/2 (r = —1/2) mass spectrum from tidal-lobe accretion 
(where Salpeter is 7 = —2.35). The higher-mass stars develop the 
steeper 7 = —2 and 7 = —5/2 (r = —3/2) mass spectra due 
to the accretion. The solid lines demonstrate the 7 = — 2 and 
7 = —2.5 slopes. 

case P is dependent on M*(0) and we cannot derive an an- 
alytical expression for the resulting mass spectrum. In the 
case where the incoming mass distribution is from the tidal- 
lobe accretion, then from §3.1 we have an expression for the 
initial mass-radius correlation (equation [l?]) . Using this, we 
can generate a mass spectrum from accretion for a cluster 
in a n oc R~ 2 configuration. Numerical simulations of accre- 
tion in stellar dominated clusters have values of 77 = 1/3 and 
C = 3/2 (Bonnell et al. 2000). Equation (||) then becomes 



M,ocM 2 R~?. (33) 

Figure ^ shows the resultant IMF for clusters accreting from 
an initial delta function mass distribution and those with 
an initially M^ 3 ^ 2 mass spectrum that accretes according 
to equation (^). In the former case, the stars, in a n oc r~ 2 
distribution, accrete with a spatially-independent accretion 
rate (M* oc M 2 ). These clusters develop an M^ 2 mass spec- 
trum. In contrast, with the initial mass segregation, the IMF 
develops into a steeper (7 ft; —2.5) IMF as the stars accrete 
in the stellar dominated regime. This is somewhat steeper 
than the analytical MZ 2 mass spectrum due to the initial 
correlation between stellar mass and position in the cluster. 

We expect stellar dominated clusters to accrete towards 
a mass spectrum with asymptotic limits of — 2.5 < 7 < —2. 
This mass spectrum is expected for the higher mass stars as 
it is they that are located in the centre of the cluster and 
have attained a higher mass due to their previous accretion 
in the gas dominated regime. Clusters that have mass dis- 
tributions that are uncorrclated with cluster position going 
into the stellar dominated phase result in mass spectra of 
7 » —2 whereas those that are mass segregated initially will 
have slightly steeper mass spectra 7 ~ —2.5. 

It is important to note that the above derivation relies 
upon all the stars stopping their accretion at the same time, 
as the final masses are a strong function of the time (equa- 
tion This should not be too strong a limitation as in the 
stellar dominated potential the accretion timescale must be 
much longer than the local crossing time as M gas << M sta rs. 
In a rich young stellar cluster, gas removal is due to the pres- 
ence of massive stars. These stars ionise the gas such that the 
sound speed (w 10 km/s) is comparable to or greater than 
the velocity dispersion. Thus the gas should be removed on 
a timescale shorter than or comparable to the crossing time, 
and all the stars in the core should stop accreting quasi- 
simultaneously. If this is not the case, then the resulting 
mass spectrum will exhibit a local maximum corresponding 
to stars that are just outside the cleared region. 

In summary, we expect that the combination of a gas 
dominated and a stellar dominated regimes and the different 
accretion physics operating in each, results in a two power- 
law IMF. The lower-mass stars have a shallower (7 « —1.5) 
slope as their mass accumulation is dominated by tidal-lobe 
accretion wheras the higher mass stars have a steeper (—2 > 
7 < —2.5) slope as their mass accumulation is dominated by 
Bondi-Hoyle accretion in the stellar dominated core. 



4 SIMULATED CLUSTER IMFS 

In order to explore the relevance of the asymptotic limits for 
the mass spectrum, we performed a number of simulations 
of accretion onto a cluster of 1000 stars. The stars are ini- 
tially of equal mass (O.IMq) and are embedded in gas which 
comprises 91 per cent of the total cluster mass The gas is 
cold such that it contains 1000 Jeans masses, 




where p is the gas density, T is the gas temperature, R g 
is the gas constant, G is the gravitational constant, and [i 
is the mean molecular weight. The simulations were per- 
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Figure 3. Histograms of the resultant mass spectra are plotted 
(in log mass bins) for a cluster containing 1000 stars that initially 
comprise 10 per cent of the total mass. The mass spectra are pre- 
sented at three different times (dotted line, t = 0.5%; dashed line, 
t = 0.8tfj; and solid line, t = 0.96%. The solid lines denote the 
expected slopes of 7 = —3/2 for low-mass stars and 7 = —2 and 
7 = —5/2 for high-mass stars. 



formed with a hybrid N-body SPH code (Bate, Bonnell & 
Price 1995) which uses standard SPH particles to model the 
gas and sink-particles to model the stars (for more details 
see also Bonnell et al. 2000). These sink-particles interact 
only gravitationally with the rest of the cluster and by ac- 
creting gas particles that come within their sink radius. This 
sink or accretion radius was taken to be the smaller of the 
tidal-lobe radius, the Bondi-Hoyle radius or the stellar sepa- 
ration. In addition, gas particles can only be accreted if they 
are bound to the star. The simulations were performed with 
either 9000 or 90000 SPH particles. The low-resolution sim- 
ulations resulted in an excess of stars in the lowest (initial) 
mass bins as they were not able to resolve the low accretion 
rates involved. The higher-resolution simulations resolved 
the accretion onto all stars and both low and high-resolution 
simulations resulted in similar high-mass mass-spectra. 

Both the stellar and gas distribution are initially uni- 
form but collapse down towards a centrally condensed dis- 
tribution. Gas is accreted by the stars and removed from 
the simulation. As the gas is accreted preferentially near 
the centre once a power-law density profile is established 
(p oc R~ 2 ), the stars soon dominate the potential there. The 
gas must then infall from further out in the cluster. The sim- 
ulations result in mass spectra that are broadly consistent 
with our analytical expectations. The low-mass stars that 
accrete their mass during the gas dominated phase display 
shallow (7 « —3/2) mass spectra whereas the high-mass 
stars that accrete the majority of their mass during the 
stellar dominated phase display steeper (—2.5 < 7 < —2) 
mass spectra. As the clusters are initially uniform and cold, 
only the end of the simulations develop a power-law density 
profile and that just before entering the stellar dominated 
regime. Because of this, the expected low-mass mass spec- 
trum does not extend over a large mass range. 

In Figure |[ we plot the resultant mass function for one 
of the higher resolution simulations, (performed with 90000 
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Figure 4. The mass accreted during the stellar-dominated phase, 
MasssDP, is plotted against total stellar mass. High mass stars de- 
rive the vast majority of their mass from the accretion in a stellar 
dominated potential whereas the majority of low-mass stars ac- 
cumulate their mass predominantly in the gas-dominated phase. 
The continuous line represents where stars would lie if all their 
mass was accreted during the stellar-dominated phase whereas 
the dotted line represents where stars would lie if half of their 
final mass was accreted during this phase. The solid histogram 
gives the mean mass accreted in the stellar dominated potential 
as a function of final stellar mass. 



SPH particles) at three different times during the evolution. 
The mass spectrum increases in breadth with time and de- 
velops into a power-law for the high-mass stars. The low- 
mass mass spectrum is broadly consistent with tidal-lobe 
accretion in a mostly uniform cloud, although only over a 
small range in mass. Due to the uniform initial conditions, 
only a fraction of the stars experience the power-law density 
profile before entering the stellar-dominated phase. This re- 
gion would extend further in mass if the cluster was initially 
centrally condensed when accretion begins. At later times, 
the high-mass mass spectrum converges to a slope between 
7 = — 2 and 7 = —2.5 consistent with the accretion during 
the stellar dominated phase, and to the Salpeter slope of 
7 = -2.35. 



4.1 Accretion dynamics 

In order ascertain whether we are correct in our interpreta- 
tion of a two-power law IMF resulting from accretion in gas- 
dominated and stellar-dominated regimes, we evaluated how 
much of the eventual stellar mass was added in each regime. 
Figure ^ plots the amount of mass accumulated by each star 
during the stellar-dominated phase against the final stellar 
mass. High-mass stars accumulate the majority of their mass 
during this phase whereas most of the low-mass stars only 
accrete a small fraction of their mass during this phase. The 
break between the two regimes occurs at approximately the 
same final mass (~ IMq) where the mass function displays 
a break between the two slopes. This supports the assertion 
that the two different power-laws in the IMF derive from dif- 
ferent physical regimes which affect how the stars accrete. 
The low-mass stars derive their mass from tidal-lobe accre- 



© RAS, MNRAS 000, [!]-§ 



Accretion and the IMF 7 



tion during the gas-dominated regime whereas the high-mass 
stars derive the majority of their mass from a Bondi-Hoyle 
type accretion that occurs in the inner parts of the cluster 
where the potential is dominated by the stars themselves. 



5 DISCUSSION 

In addition to producing the two power-law mass spectrum, 
competitive accretion naturally results in a certain degree 
of mass segregation. This arises due to the accretion in the 
gas-dominated phase where there is a strong correlation be- 
tween accretion rate, and thus the final mass, and position 
in the cluster (see equations (p^) and ( po[ ) . This direct corre- 
lation between the final mass and position in the cluster ne- 
glects variations in the initial masses and the relative move- 
ments of the stars due to their interactions. If the cluster is 
mass segregated entering the stellar dominated phase, then 
the M, oc Ml implies that the mass segregation will per- 
sist. Simulations of accretion in clusters show that the mass 
segregation does result but that there is not a one to one 
correlation between mass and radius (Bonnell et al. 2000). 
In fact, low-mass stars are located throughout the cluster, 
including in the core, but the high-mass stars are predom- 
inantly located in the central regions as is found in young 
stellar clusters such as the ONC (Hillenbrand 1997). 

It is also worth noting that although the models pre- 
sented here are meant to consider accretion onto young stars, 
they are equally appropriate for the growth of clumps in 
a molecular cloud. As the clumps evolve towards gravita- 
tional instability, they will accrete from the surrounding gas 
and this accretion will be governed by the physics described 
here. Thus, for example, the clump mass-function found by 
Motte, Andre & Neri (1998) for the p Oph molecular cloud 
could be due to the accretion by the pre-stellar clumps as 
the whole system collapses down to form a cluster. The ob- 
served 7 = —3/2 slope would imply that the whole system 
is in a p cx R~ 2 density configuration and is subvirial (dom- 
inated by the diffuse gas not in the clumps). The steeper 
slope found by Motte et al. (1998) at the high- mass end of 
the mass spectrum can be interpreted as arising from a re- 
gion of near-uniform gas density. A test of such a possibility 
is to estimate the degree of mass segregation of the clumps 
in this pre-stellar cluster system. 

Finally, it is possible that the mass spectrum for massive 
stars, M* > IOM0, is significantly different than that for 
lower mass stars if they do not form in a similar fashion. If 
massive stars cannot accrete above IOMq due to the effect 
of radiation pressure on the infalling dust (Yorke & Kriigel 
1977; Yorke 1993) but form through a merger process in 
a dense core (Bonnell et al. 1998) then the expected mass 
spectrum could be significantly different from that presented 
here. 



6 CONCLUSIONS 

Competitive accretion in young stellar clusters results in a 
two power-law mass spectrum with a slope of 7 ~ —3/2 for 
low-mass stars and a steeper slope of —2 > 7> — 2.5 for high- 
mass stars. The different slopes are due to whether the gas 



or the stars dominate the cluster potential. When gas domi- 
nates the cluster potential, the accretion rates are given by a 
tidal-lobe accretion radius and results in a dN cx MZ 3/2 dM r 
mass spectrum in an isothermal sphere distribution. The 
shallower density profile in the central regions implies a 
steeper mass spectrum, 3/2 > 7 > —4. Once the stars have 
accreted enough to dominate the potential, which occurs 
from the inside out, the appropriate accretion radius is the 
Bondi-Hoyle radius. When this happens the strong depen- 
dence of the accretion rate on stellar mass M* cx Ml results 
in a mass spectrum of dN cx M* 2 dM* to dN cx Mr 2 ' 5 dM*. 

Simulations of accretion in clusters containing 1000 
stars result in mass spectra broadly consistent with these 
analytic estimates. The high-mass stars are confirmed to ac- 
cumulate the majority of their mass in a stellar-dominated 
potential. Lastly, competitive accretion in clusters naturally 
results in mass segregation as the accretion rates are higher 
near the cluster centre during the gas dominated phase. 
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